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The weak decays of A;, ^ A + 7 and A^ ^ A + are investigated in the Standard 
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Model using hght-cone sum rules approach. The higher twist distribution amplitudes of 
A baryon to the leading confornial spin are included in the sum rules for transition form 
factors. Our results indicate that the higher twist distribution amplitudes almost have no 
' influences on the transition form factors retaining the heavy quark spin symmetry, while such 

corrections can result in significant impacts on the form factors breaking the heavy quark 
spin symmetry. Two phenomenological models (COZ and FZOZ) for the wave function of A 
baryon are also employed in the sum rules for a comparison, which can give rise to the form 
' factors approximately five times larger than that in terms of conformal expansion. Utilizing 

^ ! 

00 \ the form factors calculated in LCSR, the physical observables like decay rate, polarization 

\^ . asymmetry and forward-backward asymmetry are analyzed for the decays of Ah — > A7, 

P ■ Az+r. 

O 
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O ■ I. INTRODUCTION 

> ' 

^ ' Generally, new physics can be accessible through rare decays, where the contributions from the 

H : 

Standard Model (SM) are suppressed enough. Hence, such decays can provide an ideal platform 
to test the SM precisely as well as to bound new physics parameters stringently. Rare decays 
involving b ^ s flavor changing neutral current (FCNC), which are forbidden at the tree level in 



the standard model, can only be induced by Glashow-Iliopoa 



£s-Maiani mechanism |l| via loop 



diagrams. The Cabibbo-Kobayashi-Maskawa (CKM) matrix ^, |3|] elements can be determined 
quantitatively from 6 — > s rare decays, — mixing [4] etc., which will test its unitarity under 
the requirement of the SM. 

It is well known that the inclusive decays are relatively robust theoretically, since the decay rate 
can be systematically and reasonably approximated by the decay of a free b quark into light quarks, 
gluons and photons; but the counterpart of experimental measurements are quite difficult. On the 
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contrary, the experimental investigations of exclusive decays are comparably easier; while the the- 
oretical analysis of these modes encounter formidable difficulties due to lack of good understanding 
of QCD at low energy. As a matter of fact, there have been extensive studies on the exclusive 
decays oiB^Ml+l- Q, 13, Q, 13, fl Q] , 
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13, 
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B y(A)7 

the literature with varying degrees of theoretical rigor and emphasis. Unlike mesonic decays, the 
investigations of FCNC b ^ s transition for bottom baryonic decays Af, — > A7 and A;, — > Al~^l~ are 
much behind because more degrees of freedom are involved in the bound state of baryon system 
at the quark level. It should be pointed out that such baryonic decays can offer the unique ground 
to extract the helicity structure of effective Hamiltonian for FCNC b ^ s transition in the SM 
and beyond, which is lost in the hadronization of meson case. From the viewpoint of experiment, 
the only drawback of bottom baryon decays is that the production rate of A^ baryon in b quark 
hadronization is about four times less than that of the B meson. 

The polarization asymmetry and forward-backward asymmetry in rare decays are very sensitive 
to the new physics effects beyond the SM 
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2^. For instance, the 



polarization asymmetry of A baryon in A^ ^ A + decays is dependent heavily on the right- 
handed current j^], which is suppressed in the SM. However, before claiming the new physics signal, 
one should have a serious and comprehensive scrutiny of discrepancies between experimental data 
and theoretical predictions within the SM to confirm whether they are indeed due to contributions 
from new physics or missing effects in the SM, such as higher power and next-to-leading order 
corrections, soft gluon effects and so on. 

Currently, there has been some studies in the literature on A;, ^ A transition form factors 
ranging from phenomenological models including the pole model (PM) [30], covariant oscillator 
quark model (C0QM)[3li, MIT bag model (BM)[33] and non-relativistic quark model [331] to QCD 
sum rule approach (QCDSR) sj] and also perturbative QCD (PQCD) approach 35]. It can be 
observed that the available theoretical predictions vary from each other and can differ even by the 
orders of magnitude. For instance, the decay width of A^ A7 in PQCD approach is about one 
to two orders smaller than that obtained within other frameworks. 

It is suggested that the soft non-perturbative contribution to the transition form factor 
can be calculated quantitatively in the framework of light-cone sum rules (LCSR) approach 
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, which is a fully relativistic approach and well rooted in quantum field theory, 
in a systematic and almost model-independent way. As a marriage of standard QCDSR tech- 



nique 
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43| and theory of hard exclusive process 
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5l|, LCSR 



cure the problem of QCDSR applying to the large momentum transfer by performing the op- 
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erator product expansion (OPE) in terms of twist of the revelent operators rather than their 



dimension 



52] • Therefore, the principal discrepancy between QCDSR and LCSR consists in that 



non-perturbative vacuum condensates representing the long-distance quark and gluon interactions 
in the short-distance expansion are substituted by the light cone distribution amplitudes (LCDAs) 
describing the distribution of longitudinal momentum carried by the valence quarks of hadronic 
bound system in the expansion of transverse-distance between partons in the infinite momentum 
frame. Phenomenologically, LCSR has been widely applied to investigate the transition of mesons 



and baryons in the recent years, such as bottom meson decays 
netic form factors of nucleon 
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of Ab plvi [63], Ac ^ Mvi 
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53|, 



6l[ |. pion electroproduction 
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57|, 



5a], electromag- 



62], semi-leptonic weak decays 



641 ]. and can make consistent predictions with experimental data 
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N transition 

n 



661 ]. radiative 



successfully. The generalization to the 7*iV ^ A transition 

I 

decays of decuplet to octet baryons |67|] and meson-octet baryon couplings J68|] have also been 
worked out. 

LCDAs are the fundamental non-perturbative functions in the LCSR, which describes the 
hadronic structure in rare parton configurations with a fixed number of Fock components at small 
transverse separation in the infinite momentum frame. There have been continuous interests con- 
centrating on the research of pre-asymptotic corrections to the distribution amplitudes of hadrons 
in the exclusive reactions over two decades, in an attempt to accommodate the experimental data, 
such as the electromagnetic form factors of pion and nucleon 
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par ticular, the LCDAs of 
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mesons including higher twist have been investigated extensively 
contrast to the meticulous and intensive study of mesonic LCDAs, the corresponding task for the 
baryonic counterpart received less attention in the past due to the complexity of the inner struc- 
tures for baryon. Fortunately, an exciting attempt to study the LCDAs of baryon systematically 
was carried out in [3] for the first time, where all the transverse degrees of freedom are integrated 
out and the higher twist distribution amplitudes contributed by the light-ray operators involving 
"minus" components of the quark field operators were investigated in detail. Subsequently, an 
alternative approach to classify the three-quark light-cone amplitudes of proton was put forward 
in 



771, 



79f ]. where the transverse degree of freedom of partons is retained and the minus com- 



ponent of field operators are eliminated in favor of the plus and transverse ones with the help of 



equations of motion. It is mentioned in Ref. [8(1] that both techniques are probably equivalent. In 
the present work, we would like to follow the prescription presented in 13] for LCSR approach to 
study the rare decay of A^ ^ A-I-7 and A^ — > A + which can be also generalized to the decays 
of Af, to heavier A baryons 8l|, [82] , once the distribution amplitudes of corresponding baryons are 
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available. 

The organization of this paper is as following: In section II, we collect the effective Hamiltonian 
responsible for FCNC transition of 6 — > 57 and b — > sl^l~ . Parameterizations of various hadronic 
matrix elements {A\sTib\Ai,) with Fj being all the possible Lorentz structures are also presented 
here. Section III contains the derivation of LCSR for both A;, — > A + 7 and A^ — > A + decays 
up to the leading order of Os with the corrections from higher twists. We find that the relations 
between various transition form factors existing in the heavy quark effective theory (HQET) and 
large energy effective theory (LEET) are well respected within the framework of LCSR approach. 
The numerical analysis of LCSR for transition form factors are displayed in section IV, where 
the corrections from the distribution amplitudes of higher twist and comparisons with the results 
that derived in other methods are investigated at length. In section V, we apply the available 
form factors to study the invariant mass distribution of lepton pair, A polarization asymmetry and 
forward-backward asymmetry of A;, ^ A-|- l^l^ together with the decay width and A polarization 
asymmetry of A^ — > A7. We confirm the finding in Ref. {34, 



83( 1 that the A polarization asymmetry 



of A;, A7 is completely independent on the hadronic transition form factors, which makes it 
a good quantity to study the chiral structure of 6 — > s effective Hamiltonian and probe the new 
physics beyond the SM without receiving any pollution from strong interactions. Numerical results 
of such observables and comparisons with that obtained in other theoretical methods, in particular 
QCDSR and PQCD approach, are also included in this section. Section VI is devoted to our 
conclusions. 



II. EFFECTIVE HAMILTONIAN AND PARAMETERIZATIONS OF MATRIX 

ELEMENT 

A. Effective Hamiltonian for 6 — > s transition 



Integrating out the particles including top quark, and Z bosons above scale // = 0{mi)) , 



we arrive at the effective Hamiltonian responsible for the b sl^l transition 

6 



84| 



H,ff{b^sl+l-) 



Gp 
2y/2 



i=l 



(1) 



where we have neglected the terms proportional to 
the complete list of the operators can be found in 



:uhV*g on account of \ VubV*g/VtbV*g\ < 0.02 and 
84( 1 . In the SM, Wilson coefficients Cj at scale 
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TABLE I: Numerical values of Wilson coefficients in the NDR scheme at scale /i = rrib 



Ci 


C2 


C3 


C4 


C5 


Ce 


C77 


CsG 


C9 


C'lo 


-0.218 


1.107 


0.011 


-0.026 


0.007 


-0.031 


-0.305 


-0.15 


4.344 


-4.669 



/i = nil, calculated in the naive dimensional regularization (NDR) scheme are collected in Table HI 
In terms of the Hamiltonian in Eq. ([1]), we can derive the free quark decay amplitude for 
b — > sl^l^ process as 



Gf 



2i 



2^/2 vr 



(2) 



We should emphasize that the Wilson coefficient Cio does not renormalize under QCD corrections 
and hence is independent on the energy scale ^ ~ 0(mf,), since the operator Oio can not be 
induced by the insertion of four-quark operators due to the absence of Z boson in the effective 
theory. Moreover, the above quark decay amplitude can also receive additional contributions from 
the matrix element of four-quark operators, X]^=i(^^^~s|Oi|6), which are usually absorbed into 

n) . To be more specific, we can decompose Cg^^ (fi) into the 



the effective Wilson coefficient Cq^^ 



following three parts 
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9l| 



Cl"{fl) = Cgifl) + Ysd{z, s') + Yld{z, s'), 



(3) 



where the parameters z and s' are defined as z = rric/mb, s' = q^/m^. Ysd{z-,s') describes the 
short-distance contributions from four-quark operators far away from the cc resonance regions, 
which can be calculated reliably in perturbative theory. The long-distance contributions IldI-^, s') 
from four-quark operators near the cc resonance cannot be calculated from first principles of QCD 
and are usually parameterized in the form of a phenomenological Breit-Wigner formula making 



use of the vacuum saturation approximation and guar 
for Ysd{z-,s') and Ylo{z-,s') can be written as 



23 



c-hadron duality. The manifest expressions 



24|, 



Ysd{z,s') = h{z,s'){?,Cl{^l) + C2{^l) + ^c■i{|I) + c^{^l) + ?,c^{^Ji) + Cf,{^J)) 

-^/i(0, /)(C3(/x) + 3C74(^)) + ^(3C3(/u) + C74(^) + 3C5(^) + ^(m)), (4) 



6 



YLDiz,s') 



ot 



-(3Ci(/i) + C2(/i) + 3C3(m) + C4(/i) + 3C5(^) + C6(m)) 



j=ip,ip' 



7rr{j l+l-)Mj 



(5) 



with 



/i(0,s') 



--Inz + — + -X - -(2 + x)|l - xl^/^ ^ 
9 27 9 9^ ^' ' ^ 

8 8 mfe 4 4. 

— m ms H — ZTT . 

27 9 /i 9 9 



In 



2 arctan 



ivr for x = Az^ / s' < 1 
for X = Az^ / s' > 1 



(6) 



Here Mj{Tj'*) are the masses (widths) of the intermediate resonant states and T{j l^l ) denote 
the partial decay width for the transition of vector charmonium state to massless lepton pair, which 



can be expressed in terms of the decay constant of charmonium through the relation 



(7) 



The phenomenological parameter kj in Eq.(l5|) is introduced to account for inadequacies of the 
factorization approximation, which can be determined from 



BR{Ai, KJ/tl) M+r) = BR{Ah AJ/^) ■ BR{J/^P l+T). 



(8) 



The function Wj(g^) introduced in Eq.(l5|) is to compensate the naive treatment of long distance 
contributions due to the charm quark loop in the spirit of quark-hadron duality, which can over- 
estimate the genuine effect of the charm quark at small remarkably ^. The quantity u}j{q^) 
can be normalized to ujj{M^ ) = 1, but its exact form is unknown at present. It has been shown 
that there is a suppression for ujj{q^) when going down to = with w(0) < 0.13 (and could 
be even smaller) [98t], hence we can neglect the resonant contributions for b ^ safely. Since 
the dominant contribution of the resonances is in the vicinity of the intermediate ipi masses, we 
will simply use u]j{q^) = 1 in our numerical calculations. Moreover, the non-factorizable effects 



94 



99| from the charm loop can bring about further corrections to the radiative b ^ S'y 
transition, which can be absorbed into the effective Wilson coefficient C^-^-^ as usual. Specifically, 



the Wilson coefficient C^-^^ is given by 



22 



(9) 



^ For a more detailed discussion on long-distance and short-distance contributions from the charm loop, one can 
refer to references ^ HI, HI, [H, [H, [oa, [13 ■ 
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with 



^^/^^(Gi(xt) - 0.1687) - 0.03C2(^) 



Gi(x) 



8(x - 1)3 



where ry = as(mvy)/as(/^), = ml/ni^r, C'^_ 



>S7 



4(x-l)4' 
is the absorptive part for the h 



sec 



(10) 
(11) 



rescattering and we have dropped out the tiny contributions proportion to CKM sector VubVus- 
Similarly, the free quark decay amplitude for 6 — > 57 can be written as 84i |: 



Gp 



M{b^sj) = ^VtbVZ^C'/f{fi)sa^,q''{mbR + msL)bF 
We stress again that one should add a term 



(12) 



[3C7i(/i) + C72(/i) + SC-iifi) + C4(/i) + 3C5(/u) + Ceif,)] — ^ u;,(0)A: 



7rr(j ^ l+l-)Mj 



j=ip,ip' 



tot 



(13) 



to the effective Wilson coefficient C^^^ , if the long-distance contributions from the charm quark 
loop in the resonance regions are included. As mentioned above, this type of effects are suppressed 
heavily both by the Breit-Wigner factor ~ Fj/Mj and phenomenological parameter w(0), and hence 
are neglected in the later discussions. 



B. Parameterizations of hadronic matrix element 



With the free quark decay amplitude available, we can proceed to calculate the decay amplitudes 
for Afe — > A7 and A;, Al^l~ at hadron level, which can be obtained by sandwiching the free quark 
amplitudes between the initial and final baryon states. Consequently, the following four hadronic 
matrix elements 



{A{P)\sj^b\Ab{P + q)) , {A{P)\sj^-f5b\Ab{P + q)), 
{A{P)\saf,,b\Ab{P + q)) , (A(P)|sa^,756|Ab(P + g)), 



(14) 
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need to be computed as can be observed from Eqs. ([2]) and ([IS]). Generally, the above four matrix 



elements can be parameterized in terms of a series of form factors as 



22 



26 



271, 



28 



29| 



(A(P)|s7^6|A,(P + q)) = A(P)(<7i7^ + ffsia^.g" + 93q^)MP + q), (15) 

(A(P)|s7^756|A,(P + q)) = A(P)(Gi7^ + G2ia^uq'' + G^q^)^MP + <l), (16) 
{A(P)\sa^^b\Ah{P + q)) = A(P)[/iicj^^ - i/i2(7/.9i. - 7^9^) 

-ihzilt^Pu - luPt.) - ih^iPf^qu - Puqi.)]Ai,{P + q), (17) 
(A(P)|scj^^756|A;,(P + q)) = A{P)[Hiaf,^ - iH2{-f^q^. - -iuq^^ 

-iHsi^f^Pu - luPf.) - iH^P^^q, - Puqf,)h5AbiP + q), (18) 

where all the form factors gi, Gi, hi and Hi are functions of the square of momentum transfer 
q'^. Rewriting the Eqs. (I17lll8p . we have the baryon matrix elements for the dipole operators as 
following 

(A(P)I w^,(7^6|Afe(P + q)) = MP){fa^ + f2ia^.q'' + hq^)Ab{P + q), (19) 
(A(P)I w^,75(7"6|Afe(P + q)) = A(P)(Pi7^ + F^ia^.q" + Fsq^)^MP + q), (20) 

with 

_ 2/i2 - /t3 + h^^rriA^ + ttt-a) ^2 ^2i) 
j.^ _ 2hi + /i3(mA - TWA J + h4q'^ ^^^^ 

h = ^^^^^h, (23) 

q 

P 2^2 - ^3 + Hi{mA, - mA) „2 /^.^ 
-^1 = 2 ^ ' 

2Hi + Hsimx + ^aJ + H^q'^ 
^2 = 2 ' ^^^^ 

P3 = !!^^^P, (26) 

q^ 

It should be emphasized that the form factors and P3 do not contribute to the decay amplitude 
of A;, ^ A + l^l^ due to the conservation of vector current, namely q^lj^l = 0. Concentrating 
on the radiative decay of A;, A7, we then observe that the matrix element of magnetic penguin 
operators can be simplified as 

{A{P)\sia^,q''b\Ah{P + q)) = f2{0)A{P)ia,,,q'' Ak{P + q), (27) 
{A(P)\sia^,^5q''b\Ah{P + q)) = F2{Q)A{P)ia ^.^^q" Ab{P + q) . (28) 
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For the completeness, we also present the parameterizations of matrix elements involving the 
scalar sb and pseudo-scalar S756 currents, which can be obtained from the Eqs. (llSp and ()16p by 
contracting both sides to the four-momentum 



{AiP)\sb\A,iP + q)) 
{A{P)\s-f,b\A,{P + q)) 



1 ^ 



A{P)[gi{mA, - ruA) + g3q^]Ab{P + q), 



(29) 



rub + nis 

^ -A{P)[Gi{mA,+mA)-G3q^]jMP + Q)- (30) 



rrib - TUs 

The two independent form factors ^1 and ^2 in HQET are defined as 

{A{P)\bTs\Ab{P + q)) = A(P)[6(g')+ ^6(g')]rA5(P + 



(31) 



with r being an arbitrary Lorentz structure and being the four-ve 
Eqs. (fTMTe]) . (fT9ll20]) and the Eq. one can easily find 

^2 



22 



26 



ocity o f Ap baryon. Comparing 



271, 



28 



fi 

h 
F3 



-6 



ruA 



F2 = gi = Gi=ii + 6, 

ruA - ruAf, , 



mA + ^ 



-6 



52 — G2 — §3 — G3 



6 



(32) 
(33) 
(34) 
(35) 
(36) 



It is known that Eq. (jSip is successful at zero recoil region (with large g^) in the heavy quark 
limit. As for the large recoil region, the large energy effective theory implies that the form factors 
are independent of both energy of light hadron (E) and the heavy quark mass (nib) with the 
assumption of Feynman mechanism, which indic ates that Eq. ()3ip is still well defined owing to the 
tiny effects from l/m^, 1/E and corrections 100(]. 
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III. LIGHT-CONE SUM RULES FOR THE TRANSITION FORM FACTORS 

The most general decomposition of matrix element for three-quark operator with lig ht-like 
separations x'^ ^ between vacuum and A baryon state at tree level can be written as 7a] 

4(0|e*^''=<(aix)4(a2x)4(a3x)|A(P)) 

= SiuiaC^p (75A)^ + S2mlC^p {^75^)^ + Vimf, (tsC)^^ + Viml (75^)^^ (^A)^ 

+Vi {pC)^p (75A)^ + VsmA {pC)^p {^^,k\ + V^m^ (7^C)„^ (t'^TsA)^ 

+V4mi {^C)^p (75A)^ + Vsmi {^^C)^^ (ia'^'^x^^.A)^ + V,ml (/C)„^ (;^f75A)^ 

+^1 75C)„;3 A^ + A2mA {Pl^C)^p + ^s^A (7m75C)^^ (7^A)^ 

+Aiml {^75C)^p A^ + A^ml (7/.75C)„/3 {io^^'x^K)^ + Ami {^-i^C)^^ {^K)^ 

+Ti {P''ia^,C)^p (t'^TsA)^ + TsmA {x^^ ia ^,C) (75A)^ + TgrnA (^r^.C)^^ (^r^^TsA)^ 

+r4mA (P'^a^.C)^^ (a^^Xp75A)^ + Tsmi (x'^ia^.C)^^ (t^TsA)^ 

+%ml (x^P'^ia^.C)^^ (^75A)^ + Tymi {a^.C)^^ {a^'' ;6j,A)^ 

+%ml {x'a^,C)^p (a^''xp75A)^ , (37) 

in view of the Lorentz covariance, spin and parity of A baryon. Here, each of the 24 invariant 
functions J-i, {T = S,V,V,A,T) is a function of the scalar product P ■ x and the parameters Oj 
denote coordinates of valence quarks. 

The "calligraphic" functions J^i defined above do not have a definite twist. Decomposing the 
quark field operators into "plus" and "minus" components j70] and utilizing the equation of motion, 
the invariant function J^i can be expressed by the distribution amplitudes \q4]: 

.Ai = ^1 , 2p- xAi = -Ax + ^2 - A , 
2Az = A-i, 4p- XA4 = -2Ai -A3-A4 + 2A5 , 
4p ■ XA5 = A3 - ^4 , {2p ■ x fAe = Ai-A2 + A3 + A4-A5 + Ae, (38) 

where the twists of distribution amplitudes Ai{i = 1,...,6) can be easily determined by power 
countering as collected in Table [III Only the the relations associated with axial- vector distribution 
amplitudes Ai of A baryon are involved in the sum rules of form factors responsible for A;, ^ A 
transition, since the spectator quarks [ud] in the Af, baryon maintain their position at the origin 
up to the accuracy of leading order in Og ■ 

Generally, the distribution amplitudes depend on the renormalization scale and can be expanded 
in contributions of light-ray operators with increasing conformal spin, where the constraints on 
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TABLE II: Twist classification [76| of the distribution amplitudes given in Eq. (|38p. 





twist-3 


twist-4 


twist-5 


twist-6 


axial- vector 


Ai 


A2 , A3 


Ai , A5 


Ae 



oper ator mixing and equation of motion owing to conformal symmetry of the QCD Lagrangian 



10l[ | are also taken into account. To the leading logarithmic accuracy, the different conformal 
partial waves do not mix with each other, since the renormalization group (RG) are driven by 
tree-level counter terms and the conformal symmetry is well respected. Besides, it is known that 
QCDSR tends to overestimate the matrix elements corresponding to the higher conformal spin 
operators considerably from the experience of both calculating the pion electromagnetic form factor 

102 . 103 1 and nucleon form factor [sol. Therefore, following the Ref. we would like to adopt 
the distribution amplitudes of A baryon to the leading conformal spin accuracy in the present work. 
The explicit forms of distribution amplitudes for A baryon to the accuracy of leading conformal 



spin can be expressed as 



64| 



with 



Ai{xi,X2,X3) = -120xiX2X3(/)3, 

^2(2:1, 3^2,3:3) = -24X1X204, 

yl3(xi,X2,X3) = -12X3(1 - X3)^4, 

^4(xi,X2,X3) = -3(1-X3)0^, 

yl5(xi,X2,X3) = -6X3(/>^, 

Aeixi,X2,xs) = -24, (39) 

4 = 4 = -h, <Pl = 4 = -^^, < = i^l = ^^. (40) 

The estimations of non-perturbative parameters /a and Ai in the framework of QCDSR approach 
have been presented in Ref. ^] and to make the paper self contained, these are collected in the 
Appendix [XI 

With the LCDAs of A baryon available, we are now in a position to derive the sum rules of 
transition form factors which are responsible for A^, A7 and A^ M~^l~ decays. The basic 
object in LCSR approach is the correlation function where one of the hadron is represented by 
the interpolating current with proper quantum number, such as spin, isospin, (charge) parity and 
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so on; and the other is described by its vector state manifestly. Information on the hadronic 
transition form factor can be extracted by matching the Green function calculated in two different 
representations, i.e., phenomenological and theoretical forms, with the help of dispersion relation 
under the assumption of quark-hadron duality. 



A. Light-cone sum rules for Af, ^ A + 7 



We consider the correlation function associating with — > A + 7 decay determined by the 
matrix element 



/ 



z^T,{P,q)=iz^ / d4xe-^''-(0|r{jA,(0)j.(x)}|A(P)), 



(41) 



between the vacuum and A baryon state |A(P)). Generai 



IQJ, 



ina 



l y, choice of i r iterp olating current for 



106 



ml 



1081 . 



109l |. The interpolating 



baryon, in particular the light baryon, is not unique 
field jAjj(O) denoting the A;, baryon and the transition current ju{x) describing the transition of 
Afe — > A7 are given by 

e^^'=K(0)C75 /d^(0)] ^h\Q) , 



jy{x) = z6(x)cr^^(l - 75)g^s(x) . 



(42) 



In addition, the vacuum-to-baryon matrix element for the interpolating current can be parameter- 
ized as follows 



(oiiA,(o)iA,(p')) = /A,(^-n 



(43) 



Inserting the complete set of states between the currents in Eq. (j4ip with the same quantum 
numbers as A;,, we can arrive at the hadronic representation of the correlator 



(0|jAjAfe(PO)(A,(PO|z- j.|A(P)) ^ ^ {Q\j^MP')){HP')\z^jmP)) 



mi -P' 



ml - P" 



(44) 



where P' = P + q and we have separated the contributions from the ground state and higher states 
corresponding to the A;, baryon channel. Combining the Eq. ()19l I20p . ()431 1441) and summing over 
the polarization of baryon, we can obtain the phenomenological representations of correlation 
function as 

{z-P')^ r 



z'T, = 2/a, 



'ml -P 



•i2 



fl ^+f2^4-Fi ^75 - F2 



A(P) + ..., 



(45) 



where the ellipses stand for the terms proportional to the higher power of 1/P in the infinite 
momentum kinematics P ~ 00, q ~ const, z ^ 1/P and the contributions from the higher states 
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FIG. 1: The tree level contribution to the correlation function (HT]) . where the thick solid line represents the 
heavy b quark. 



of A;, channel. The advantage of contracting the correlation function with z'^ is that the Lorentz 
structures can be simplified owing to the null contributions from on the light-cone. 

On the theoretical side, the correlation function (j4ip can also be computed in the perturbation 
theory with the help of OPE technique at the deep Euclidean region: P'\ g2 ^ _q2 ^ 0. To the 
leading order of as, the correlation function can be calculated by contracting the bottom quark 
field in Eq. (j4ip and inserting the free b quark propagator 



shown in Fig. [TJ 



lie 



llll | and can 



The full quark propagator also receives corrections from the background field 
be written as 

+ / ,, ^x^G^'^M7,], (47) 

where the first term is the free-quark propagator and G-^^ = C^^^T^ with Tr [T°T^] = ^(5"*. Inserting 
the second term proportional to the gluon field strength into the correlation function can result in 
the distribution amplitudes corresponding to the higher Fock states of A baryon. It is expected 
that such corrections associating with the LCD As of higher Fock states do not play any significant 
role in the sum rules for transition form factors 112l |. and therefore can be neglected safely in the 



presented work. 

Substituting Eq. (|37p into Eq. (j46p and performing the integral in the coordinate space, we can 
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achieve the correlation function in the momentum representation at the quark level as 

.'T„ ^-2(..P)v^(i-,.)A(p){/^. (j-;';;:l„. 

+M / ax^Xn j^-r 

+2M^ /'dx ^3 -^1(^3) - Mx2,) + islrcs) + i4(x3) - A^jx^) + ielxs) ] 
io ^ ^ [(x3P + g)2-m2]3 / 



+2q^{z-Pf /(l + 75)A(P)|M^'dx3X3- 



-il(x3)+i2(x3)-i3(x3) 
[(x3P + g)2 -m2]2 



^.M^ [\ 2 MXZ) - Mx^) + ^(3^3) + M^^) - MX^) + -46(^3) 1 , i.^. 

Jo [(x3P + g)^ -m^]3 J 

where subleading terms in the infinite momentum frame kinematics denoted by the ellipses are not 
included. The distribution amplitudes with tildes are defined as 

F(x3) = / dxg / (ixiF(xi, 1 — xi — X3, X3) , 
Ji Jo 

F(x3) = / (ixg / dxg / dxiF(xi, 1 - xi - X3,X3), (49) 
Ji Ji Jo 

originating from the partial integral in the variable X3 in order to eliminate the factor 1 /P ■ x due 
to the insertion of distribution amplitudes in Eq. ([381) . 

For the convenience of matching the correlation in QCD representation and hadronic level, the 
Eq. ()48p is usually written in a form of dispersion integral as 



where is defined as = —q^. With the assumption of quark-hadron duality, the higher states 
in the Ab channel can be given by the same dispersion integral only with the lower bound replaced 
by the effective threshold parameter sq. Besides, the Borel transformation is commonly introduced 
in the standard procedure of sum rules approach for the sake of compensating the deficiency due to 
the approximation of quark-hadron duality. Finally, we can obtain the sum rules for the transition 
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form factors of the Af, — > A7 decay as 



s/M 



Ml 



il(x3) +12(^3) -i3(x3) 



+ 



xoe 



-so /Ml 



M2 



Ml 



-M^xo-^\ 



Ai{xo) - A2{xo) + isCxo) + Mxq) - A^{xo) + A^{xo) 
xq ( Ai(xo) - ^2(3:^0) + A^{xq) + Ai{xQ) - A^{xq) + Afi{xQ) 







and 

fAj2{q')e 



-mlJMl 



M2 

Ml 
M^ 
'Ml 



M^xle-'o/M': 



ml — + XqM"^ 
M2 



dxiAi{xi, 1-Xi- X3,X3) 

- 2Ai{xs) + 12(2:3) - 13(2:3) - MX3) + A^ixs) 

MX3) - Mx3) + Mxs) + MX3) - MX3) + Mxs) 

- 2Ai{xo) + A2{xo) - A3{xo) - ^(xo) + A^ixo) 



Ml 



Ai{xo) - A2ixo) + ^3(^0) + ^4(2:0) - ^5(2:0) + Ae{xo) (52) 



-M^ 



dxn 



xl[ Ai{xo) - A2{xo) + A3{xo) + A^ixo) - A^ixo) + Afi{xo) 
m^ - q"^ + XqM^ 



with 



xo 



(Q2 + _ ^2)2 + 4M2(Q2 + ^2) _ (g2 ^ _ ^2) 

2M2 ■ 



(53) 



It can be easily observed that the sum rules for form factors Fi{q^) and F2{q^) are the same 
as that for fi{q^) and f2{q^) respectively, namely, Fi(g^) = /i(g^), -^2(9^) = f2{q^)- In addition, 
we indeed find that form factors f2{<f') and F2{q^) are exactly equal to zero at zero momentum 
transfer, which indicates that they do not play any role to the radiative decay — > A + 7. 

Now we are going to investigate the asymptotic behavior of form factors fi{q^) and f2{Q^) in the 
limit oi Q"^ ^ 00. The form factors are dominated by the configurations with the recoiled s quark 
taking most of the momentum of A baryon. Employing the asymptotic distribution amplitudes 
of A baryon and expanding the sum rules given in Eqs. ([5T]) and ([52]) in l/Q^, we can get the 
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following expressions of form factors in the limit of large momentum transfer 



2, _ 1 /a M 1 



+ 



6 /a, Ml Q6 



Ml' 



ds s e 



3„,f2 -(so-»n2 )/Af| 



+ 



ds s e 



from which one can observe that f2{q^) is suppressed by l/Q^ compared to the expected asymp- 
totic beh avior f2{Q'^) ~ in terms of the analysis on the electromagnetic pion form factor 
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UJ]. It is mentioned in 59[] that the true asymptotic behavior of f2{q^) ~ can be 



reproduced as a part of O(a^) corrections, corresponding to the hard scattering mechanism, in the 
framework of LCSR approach. We also note that the pure contributions from the leading twist 
distribution amplitudes of A baryon will give rise to the asymptotical behavior f2{Q'^) ^ l/Q^^j 
which is suppressed by one more power of l/Q^ compared to that from higher twist distribution 
amplitudes. This observation also indicates that the soft form factors of Aj, — > A are dominated 
by the distribution amplitudes with "wrong helicity " instead of the leading twist ones, which is 

59| and Af, pW 



in agreement with the analysis of soft contributions for nucleon form factors 



transition form factors 



64| 



B. Light-cone sum rules for Ab ^ A + l^l 



Different from the radiative decay A;, A7, both the matrix elements given by the tensor oper- 
ator Oj-y and vector-like currents O9, Oio sandwiching between the A^ and A states can contribute 
to the decay A;, ^ A -|- As the former matrix elements can be directly borrowed from the 

radiative decay, we only need to deal with latter one within the LCSR approach. We start with 
the following correlation function 



UP,q) = iz^ / d4xe-^''-(0|r{jA,(0)jV(x)}|A(P)), 



where the current ju{x) is given by 



jy{x) = b{x)jy{l - 75)5(2;) . 



(55) 



(56) 



We can write the phenomenological representation of the correlator at the hadronic level simply as 

(z ■ P'? 



z'T, = 2/a 



"ml -P' 



/2 



A(P) + ...., 



(57) 



where the contributions from (73 and G3 are proportional to the higher power of 1/P in the infinite 
momentum kinematics P ~ 00, q ~ const, z ~ 1/P and hence are omitted in this paper. On 
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the other hand, the correlation function at the quark level can be calculated in the framework of 
perturbative theory to the leading order of as 



z^'T, = -2{z ■ P f ^(1 - 75) A(P) <^ / -Dx 



■f 

Jo 

/' 

Jo 



X3Ai{xi) 



{X3P + g)2 - ml 



[(x3P + g)2-m2]2 



,0)1^4/ , ^AI{X3)-A2{X3)+A3{X3)+A^{X3)-A5{X3)+A6{X3)\ 

+2M / ax3a:;3 



+2(z-P)V^(l+75 



)A(P)|m^' 



dX3 X3 



-Aijxs) + A2ix3) - A3{X3) 

[(x3P + g)2-m2]2 



+2M-' / dx3X3 
/o 



+ .... (58) 



[(x3P + g)2 - m^]3 

Matching the correlation function obtained in the two different representations and performing the 
Borel transformation with respect to the variable p2^ (>an achieve the sum rules for the form 
factors gi and ^2 responsible for — > A + decay as 



r dx3e-'/^B 
Jxo 



+ 



M2 
M4 



^ da;3Ai(a;i, 1 - xi - X3,X3) 

- 2ii(x3) + A2{xs) - A3{X3) - A^ixs) + A5{X3) 



{Ai{x3) - A2{X3) + A3{X3) + A4{X3) - ^5(0:3) + ^6(^:3) 



M-Xq( 



M2 



- 2ii(xo) + A2{xq) - A3{xo) - Ai{xQ) + A^{xq) 



Ml 



Ai{xo) - A2{xq) + A3{xq) + ^4(2:0) - M{xq) + Aq{xo) (59) 



axQ 



xl ( Ai(xo) - A2(xq) + ^3(2:0) + ^4(^0) - ^5(a;o) + 
mo - + x\M'^ 



and 



,2^,-'"L/^l 



-M 



M2 



Ml 



il(x3)+ 12(2:3) -^(2:3; 



( ^i(a;3) - ^2(^3) + ^3(2:3) + Ai{x3) - A^ixs) + Aeix3) 



+ 



B 

-sa/Ml 



ml — q'^ + XqM2 
M2 



Ai{x3) + A2{X3) - A3{X3) 



Ml 



Ax{xq) - A2{xq) + ^3(^0) + Aa{xq) - A^{xq) + Aq{xq) 



(60) 



+M^Xq 



dxQ y 



Xo{ Ai{xo) - A2{xq) + A3{xo) +A4^{xo) - A5{xq) +Aq{xq) 



m^-q^ + XoM2 
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In addition, the sum rules for form factors satisfy Gi[q^) = gi{q'^), G2{q^) = 52(9^) to the accuracy 
considered in this work. These relations will be broken down with the inclusion of higher power 
corrections in l/rrih and higher order corrections in Us- Furthermore, we can also get the following 
relations between form factors /j and gi{i = 1, 2) by comparing the revelent sum rules 



which are consistent with those derived in 
symmetry and the large energy symmetry. 



/2(g')=5i(g'), 



(61) 
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lOQ l based on the analysis of heavy quark 



IV. NUMERICAL ANALYSIS OF SUM RULES FOR FORM FACTORS 

Now we are going to calculate form factors g2{Q^) and /2(^^) revelent to the — > A + 7 and 
Afe — > A + decays n umerically. Firstly, we collect the input parameters used in this paper as 



below 



3, 



115 



116 



117^ 



(62) 



Gf = 1.166 X lO-^GeV^^^ 114,1 = 41.61+J];8o x 10"^, 

\Vtb\ = 0.9991, mf, = (4.68 ± 0.03)GeV, 

mc{mc) = 1.275t|}:|]jf GeV, "i,(lGeV) = (142 ± 28)MeV, 

mA, = 5.62GeV, mA = 1.12GeV, 

/a. = 3.9t°i X 10-3GeV2, /a = 6.0l°i x IQ-^GeV^ 

Ai = -1.3l[!:^ X lO-^GeV^, sq = 39 ± IGeV^. 

It should be noted that the normalization constants of LCDAs for A;, and A baryons, namely /aj, 
/a and Ai, are all evaluated at the scale /i = IGeV. As for the choice of the threshold parameter sq, 
one should determine it by demanding the sum rules results to be relatively stable in allowed regions 
for Borel mass M^, the value of which should be around the mass square of the corresponding first 
excited states. 

With all the parameters, we can proceed to compute the numerical values of the form factors. 
In principle, the form factors g2{q^) and f2{Q^) should not depend on the Borel masse in 
a complete theory. However, as we truncate the operator product expansion up to the leading 
conformal spin of distribution amplitudes for A baryon in the leading Fock configuration and keep 
the perturbative expansion in a, to leading order, a manifest dependence of the form factors 
on the Borel parameter would emerge in practice. Therefore, one should look for a working 
"window" , where the results only mildly vary with respect to the Borel mass, so that the truncation 
is reasonable and acceptable. 
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Firstly, we concentrate on the form factors at zero momentum transfer. For the form factor 
/2(0), we require that the contributions from the higher resonances and continuum states hold the 
fraction less than 25 % in the total sum rules and the value of /2(0) does not vary drastically 
within the selected region for the Borel mass. In view of these considerations, the Borel parameter 
should not be too large in order to insure that the contributions from the higher states are 
exponentially damped as can be observed form Eqs. (fSTI) . (f52l) . (f59l) and (|60]l and the global 
quark-hadron duality is satisfactory. On the other hand, the Borel mass could not be too small 
for the validity of OPE near the light-cone for the correlation function in deep Euclidean region, 
since the contributions of higher twist distribution amplitudes amount the higher order of 1/M^ 
to the perturbative part. With the chosen threshold value sq = 39GeV^, we indeed find the Borel 
platform £ [3.0, 6.0]GeV^ which is plotted in Fig. [2l The number of /2(0) is 0.15lJ]:[!^, 
where we have combined the uncertainties from the variation of Borel parameters, fluctuation of 
threshold value, errors of b quark mass and uncertainties from the non-perturbative parameters in 
the distribution amplitudes of A baryon together. The errors on the form factor /2(0) are estimated 
within the level of 20 % as expected by the general understanding of the theoretical framework. 
Following the same method, we can continue to estimate the numerical results for the form factor 
32(0) within the selected Borel window as displayed in Fig. [2l The contributions from the excited 
resonance and continuum state are required to be less than 10 % here. It can be observed that 
c/2(0) = l.S^oi X 10-2 GeV"^ with the given Borel window M| G [3.0,6.0] GeV^. 

Next, we can further investigate the q'^ dependence of the form factors /2 and §2 based on the 
sum rules given in Eqs. ()52p and (I60|) . The OPE for the correlation function in Eqs. ()4ip and ()55p 
near the light-cone is valid only at small or intermediate squared momentum transfer, < < 
ml — 2mhAQCD, which ensure the good stability of the sum rules in Eqs. ([5T]) . ([52]) . ([59|) and ([60|) 
with the variation of M'^. The reason is that the light-cone expansion is expected to break down 
when is approaching and hence contributions from the higher twist distribution amplitudes 
increase rapidly. In phenomenology, we extend the form factors to large squared momentum 
transfer by the double-pole model 

^^(1') = 2/ f °1 4/ 4 ^ (63) 

where denotes the form factors /2 and 52- The parameters ai and 02 can be fixed by the matching 
of form factors corresponding to the small and intermediate momentum transfer calculated in the 
LCSR approach. Our results are grouped in Table IIIIl where the values of form factors in terms 
of the leading twist sum rules (twist-3) are also presented for a compassion. 
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FIG. 2: Dependence of form factors /2(0) and (?2(0) on the working Borel window M'^ £ [3.0,6.0] GeV^ 
with the chosen threshold parameter s"^ = 39GeV^. Here the bold line denotes the contributions from the 
distribution amplitudes up to twist-6, while the thin line represents the contributions from leading twist 
ones only. 



TABLE III: Numerical results for the form factors /2(0), 52(0) and parameters ai and a2 involved in the 
double-pole fit of eq. §^ for both twist-3 and twist-6 sum rules with M| e [3.0, 6.0] GeV^, sq 39±1 GeV^, 



together with results from COZ 118] and FZOZ [1191 ] models for distribution amplitudes of A baryon in the 
leading twist sum rules and results obtained in QCDSR [34]. 



parameter 



COZ FZOZ QCDSR twist-3 up to twist-6 



/2(0) 
ai 

02 



o-74i^:o6 o-87ro:o7 0.45 



14+002 



2.01^0;}^ 2Mt",:ll 0.57 2.911^^, 



^■^^-0.02 

2 94+" " 



1 Q2+014 1 41 i 



1+0.11 



-0.18 2.26 



-0.13 
-0.08 



2.31 



+0.14 
-0.10 



52(0)(10-2GeV-i) -2.4t°i -2.8t°:| 



-1.4 



-0.471^ 



ai 

02 



2.76 



+0.16 
-0.13 



2.80 



+0.16 
-0.11 



2.16 3.40 



+0.06 
06 



+0.06 
-0.05 



+0.2 

4 



2.91 



+0.12 

-0.09 



2.05t!]:?l 2.12^0:21 1.46 



2.98; 



-0.09 



2 24+° " 

^•^+-0.13 



As can be observed from Table Hill the form factor /2 is dominated by the contributions from the 
leading twist distribution amplitudes for A baryon, and the corrections owing to the higher twist 
distribution amplitudes are less than 10 %. In contrast, the higher twist distribution amphtudes 
play an important role in the form factor §2, where the numbers predicted in the twist-3 LCSR 
even differs from that of the whole sum rules in the sign. This observation is in agreement with the 
studies on semi-leptonic decay Ac — > Al~^l 6j] and nucleon form factor 59(]. Predictions on these 
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form factors can be systematically improved by including the higher conformal partial waves and 
higher order perturbative corrections. 

Apart from the general description of distribution amplitudes based on the conformal symmetry 
of QCD Lagrangian, there exist two concrete models taking into accoun t the first a few conformal 
partial waves for the A baryons, COZ model 118l |. and FZOZ model 119l |. The expressions of 



twist-3 distribution amplitudes for COZ model A^^'^{xi) and FZOZ model Af^'^^{xi) can be 
written as 

Af^^(xi,X2,X3) = -42(j)asixi,X2,X3)[0.2G{xl + xl) + 0Mxl-0.5Qx2X3-0.2Axi{x2+X3)], (64) 

Af^^^{xi,X2,X3) = -42(l)asixi,X2,X3)[0m3{xl + + 0.376x? - 0.194x2X3 - 0.207x1 (X2 + X3)], 
<?^'as(a;i,X2,X3) = 120x1x2x3. (65) 

Substituting the above distribution amplitudes for A baryon into sum rules Eqs. (I52p and (I60p and 
repeating the same procedure, we can get the parameters for oi and 02 accounting for the form 
factor /2 and (72 as grouped in Table lllli From the table, we can see that soft contributions to the 
transition form factors /2 and (72 in the COZ and FZOZ models are approximately five times larger 
than that for the LCDAs of A baryon based on the conformal spin expansion. This is similar to 
that observed in the studies of pion form factor |ll3| . nucleon form factor 59] and also Ac — > A 



transition 



The form factors are also calculated in the QCDSR [34| using the heavy quark symmetry to 
reduce the number of independent form factors. We have translated their results with the help of 
Eq. (j33ll36p . It can be observed from Table Hill that the numbers of 92(0) obtained in QCDSR 
differ from that extracted from twist-6 sum rules on the light-cone in the sign. Besides, /2(0) in 
the framework of QCDSR is about three times larger than that given by LCSR approach presented 
here. In particular, the dependence of the from factor /2 between these two methods are quite 
different. It grows slowly with the increase of squared momentum transfer in QCDSR to reach 
0.83 at the maximal momentum transfer. But it rises drastically in our approach with the increase 
of to reach 2.3 at q^ax — ("^Ab ~ "^a)^ = 20.3GeV^. Such differences between form factors will 
lead to quite different predictions on values of decay width and forward-backward asymmetry for 
semi-leptonic decay A^ — > A/"*"/". 

The hard contributions to the form factor /2(0) involving two hard gluons exchange have also 
been investigated [ssj to the leading twist distribution amplitudes for A baryon in the framework 
of PQCD approach. The values of /2(0) is computed as (1.2 ~ 1.6) x 10^^, which is about one 
order smaller than the form factor presented here in terms of the LCSR approach. As a matter of 
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fact, large soft corrections have been observed in the n qn-le ptonic charmed meson decays in the 
perturbative QCD approach based on factorization 12o!] . Therefore, it seems that the form 
factors responsible for ^ A transition may be dominated by the soft gluons exchange between 
valence quarks inside the A^ and A baryons. 

V. PHENOMENOLOGY OF A;, ^ A + 7 AND Ab ^ A + 

Utilizing the form factors derived above, we can now proceed to perform the calculations of 
decay rate, polarization asymmetry and forward-backward asymmetry. 

A. Decay width and A polarization asymmetry of A;, A + 7 

In this subsection, we present the formulae of decay rate and polarization asymmetry of A 
baryon for A;, A + 7, the latter of which can be used to analyze the helicity structures of 
effective Lagrangian at the quark level. The four-spin vector s^^ of A baryon can be defined in its 
rest frame as 

{snr.s. = (0, e), (66) 

which can be directly transformed into the rest frame of A;, baryon by the Lorentz boost 

«'^ = (^, I+^t^Pa), (67) 
rriA -Ea + rriA 

with Pa and E'a being the three-momentum and energy of A baryon. Then, the following relations 
can be directly read 

I - ^ I - ^ P 
V- s = — — 2P • s = — p • i, (68) 

1 + 2x 

where x = m\/m\^, p is a unite vector along the momentum of A baryon and v = p\^/m\^ is the 
four- velocity vector of A^. 

Making use of Eqs. (1271 ) and (1281) . the decay width of unpolarized Ah into A with a definite 

nn 

polarization vector s can be derived as [30|, 13J) 



r(A, - A7) = ^i-jyt,msm"\\^ - -'fimi + n^imm' 



2x mi — mi 



l + - 2^ l^{v-s) 

i — x^ mt + mi 



69) 



For the conven ience of comparing with the experimental data, we can rewrite the Eq. (I69p . in the 



standard form 



1151] making use of Eq. (f68]) as 



r(Ab ^ A7) = ^ro[i + a p • s] = iro[i + a' p • |], (70) 
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TABLE IV: Decay branching ratios (BR) of A(, — > A7 calculated in the LCSR approach with distribution 
amplitudes of A baryon in terms of the conformal spin expansion with only twist-3 and up to twist-6 together 
with results from COZ and FZOZ model, and results based on the form factors from pole model [30[ , covariant 



oscillator quark model 31], heavy quark effective theory 32], MIT bag model ^], non-relativistic quark 
model |33|], QCD sum rule approach 34] and perturbative QCD appr 



Model of DAs 


twist-3 


up to twist-6 COZ 


FZOZ 




BR (xlO-^) 


63+° " 








Model 


PM 


COQM HQET 


BM 


NRQM QCDSR PQCD 


BR (xl0"5) 


1.0 ~ 4.5 


0.23 0.8 - 1.5 


0.4 


0.27 3.7 ± 0.5 (0.0043 - 0.0086) 



with 



^^\Vt,m.m"\'{l - x^fiml + m^)[^(0)]^ (71) 



32m3 



and 



2x , 2x ml- ml 
" = = 1 , ^2^2 , ^2 - (72) 



1 -|- 1 + x'^ m^ + ml 



As can be observed from E g. (I72D . the polarization parameters a and a' are free of the pollution due 



to strong interactions 



831 ] and only depend on the relative strength of left- and right- handed 



couplings between quarks. Utilizing the inputs given in Eq. (|62p . we can get the polarization 
parameters as a = 0.381 ± 0.001, a' = 0.998 ± 0.001. Any distinct deviations from these values 
would indicate the new physics beyond the SM. 

Substituting the form factor /2(0) calculated in the last section into Eq. (1711) . we can achieve the 
decay rate of radiative decay A;, ^ A -|- 7 as shown in the Table IIVI Predictions on the branching 
ratio(BR) of Af, A7 can vary even by the order of magnitude adopting different models of 
distribution amplitudes for A baryon. We also collect the results of decay rate for A5 ^ A -|- 7 
computed in other approaches in Table IIVI from which we can find that different methods give 
quite different predictions. Unfortunately, only the upper bound 1.3 x 10"^ for BR of Af, A7 
decay is available in experiment at present, so we have to wait for more experimental data to 
discriminate existing models. 
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B. Decay width and dilepton distributions of A;, ^ A + l^l 

We adopt the transition form factors calculated in the sum rules up to the twist-6 distribution 
amplitudes of A baryon. The long-distance effects originated from cc resonances on the BR(A{, 
Al^l^) are also discussed . Th e differential decay width of AZ+Z~ in the rest frame of 



baryon can be written as 



1151 1 



where u = (pA = {Pi+ PA, Pi+ and pi~ are the four-momenta vectors of A, 

and /~ respectively. Maj^„»a/+«- is the decay amplitude after integrating over the angle between 
the /~ and A baryon. The upper and lower limits of u are given by 



(El + Etf - UEf -ml- JEf - m^^' 



urmn = {E*^ + E^f - U Ef - + JEf - mff; (74) 



where El. and Ef are the energies of A and / in the rest frame of lepton pair 



ml - ml 



El = ''^ - , E! = (75) 

We achieve the invariant dilepton mass distribution for Af, A + l'^l^ (l = ^,t) with and with- 
out long-distance contributions as plotted in Fig. [31 where the appearance of large enhancement at 
the end point = for Af, ^ A + /U^/x^ is due to the factor 1/g^ involved in the operator 0^-y and 
h{z,s') term in Wilson coefficient Cg'^'^^in the effective Hamiltonian for b transition. The 

dilepton mass distribution of A;, ^ A + (/ = /i, r) without long-distance contributions peaks 
at higher invariant masses in Fig. [3l which is much different from that given in [Q]. The reason 
is that the form factors grow more drastically with the increase of squared momentum transfer 
in LCSR than that obtained in QCDSR, which is also found in the semi-leptonic decay A^ — > pW 

m 



u, r) as displayed in 
Mi and PM for a 



Integrating Eq. (j73p . we get the branching fractions of A^ ^ A + l^l^ (l 
Table |Vl together with the results computed in the framework of QCDSR 

comparison. The QCDSR result {2^ for decay rate of A;, A + r+r^ without long-distance effects 
is much smaller than that presented in this work, which can be attributed to the larger transition 
form factors predicted in LCSR due to the sensitive dependence on the momentum transfer q^. 

We are now ready to investigate the effects of magnetic penguin operator Oj^ in the semi- 
leptonic decay A^ — > A + l^l^ (/ = /u, r). Without the operator Ojy, the BR(Ah — > Al^l^) is 16 
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5 10 15 20 14 16 18 20 

q^(GeV^) q^(GeV^) 



(c) (d) 

FIG. 3: The differential width for the M^l^ [l = fi,T) decays as functions of with long-distance 

contributions (a, b) and without long-distance contributions (c, d) 



% and 27 % smaller for final states being ^~ and T"'"r~, respectively. However, the branching 
fractions can be 38 % and 57 % smaller for these two modes, respectively, if we adopt the same 
magnitude of Cj^^ as that in the SM but with an opposite sign. This confirms the conclusion of 



ref. 24] that BR(Ab AZ"*"Z~) can serve as a promising quantity to explore the new physics effects 
as well as constrain parameter space of various models beyond the SM. 

Roughly speaking, the long-distance effects on the decay rate for semi-leptonic decays of — > 
A -|- can be given by 

BRLDi^b ^ A + l^r) = BR{Ab ^ A + X BRi^^t l+T) . (76) 
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TABLE V: Decay branching ratios (10^^ ) of Af, ^ A + {I = fi,T) with and without long-distance 
contributions based on the form factors calculated in this work using light-cone sum rules and that from 
QCD sum rules and pole model 



Model 


Afc ^ A + 




Ab ^ A + 






without LD with LD 


without LD with LD 


LCSR (this work) 






9 T+2.3 
^•-^-0.6 




QCDSR [24] 


2.1 


53 


0.18 


11 


PM [24] 


1.2 


36 


0.26 


9.0 



Utilizing the experimental data on the leptonic decays of J/ip and ^'(25) |ll5l ]. 

BR{J/^ = (5.93 ± 0.06)% , BR{^{2S) /U+^u") = (7.3 ± 0.8) x IQ-^ , 

Si?(^'(25) ^ r+r") = (2.8 ± 0.7) x IQ-^ , 

we can achieve the following relation 

BRld{H ^ a + /i+/x-) BR{J/^ fi+fi-) + BR{^{2S) fi+i^-) 



BRLDi^b^ A + T+T-) BR{^{2S) ^T+T-) 

= 23.8, (77) 

where the assumption BR{Ab A + J/ip) = BR{Ab ^ A -|- ^(25)) has been used in the above 
derivations. We can see that this naive estimation is consistent with that computed in this work, 
but the number got in Ref. [2^ is somewhat smaller. 

We then concentrate on the role of long-distance effects played in the decay rates for both 
final states and t^t~ within LCSR approach. As can be seen. A;, ^ A -|- ^^/U^ decay is 

dominated by the long-distance contributions in the vicinity of the cc resonance region, which is 
approximately 6 ~ 7 times larger than that from short-distance contributions. In contrast to the 
case of final states with the long-distance contributions are almost the same size as that 

from short-distance in A;, — > A + r+r" decay, which differs from that in the other two methods. 



C. A polarization asymmetry of A^ ^ A + l^l 

Similar to the radiative decay A;, ^ A -|- 7, we can write the differential decay width for 
Ab ^ A -|- with respect to the squared momentum transfer in the following form 

dT{Ab Al+l- 



A{q')+pA,-sB{q'). (78) 
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The dependence of the function describing the polarization of A baryon is plotted 

in Fig. H] for both the cases with and without LD contributions. The function A{q'^) is just the 
invariant dilepton mass distributions discussed in the previous subsection. 

Integrating the differential decay width over g^, the decay width of with polarized 

A baryon can be written in the form of Eq. (|7U|) 

r(Afe ^ A/+/-) = iro[i + a p • s] = iro[i + «'?•!], (79) 
where Tq is the total decay width of A^ ^ A + l^l^ 

/Imax 
A{q^)dq^, (80) 



28 



TABLE VI: Polarization asymmetry parameter a and a' of A baryon with and without long-distance con- 
tributions based on the form factors calculated in light-cone sum rules 









Ab^A- 


f T+T" 




without LD 


with LD 


without LD 


with LD 


a 
a' 


'J-«5'J-0.02 


-o.501°:°t 


n oa+0 03 

LI-OD_o.o3 


n 07+0.03 
'-'•^ ' -0.03 

"•'-•'^-0.03 



with q^-^ = Amf and q^ax — ("^Aj — "^-a)^- Besides, a and a' are the polarization asymmetry 
parameters, whose manifest expressions can be expressed in terms of function B{(f') as 



.2,1 ^ „„2 



2 . -Ea To y„2 mA 



with IpaI = Y ( — ^2mA )^ ~ "^A being the magnitude of three-momentum for A baryon in the 

rest frame of A;,. The numerical results of polarization variables a and a' are grouped in Table IVTl 
As can be seen, the long-distance contributions have small effects on the polarization asymmetry 
for decay of A^ ^ A -1- r+r^, however, they can modify the polarization asymmetry for the mode 
Af, ^ A -|- /i^^^ remarkably. The reason is that only cc resonances above (including) the ^(25') 
threshold can have influences on the long-distance contributions for A5 — > A -|- r+r", which are 
much suppressed compared with that from J in the decay of — > A -|- /i+^^ . The value of a 
corresponding to the channel A^ — > A -|- fi^fj,~ is consistent with — 0.54+qq4 that calculated in 
within the error bar. 



D. Forward-backward asymmetry of Ab A -I- 

For the illustration of forward-backward asymmetry, we consider the following double patrial 
differential decay rates for the decays of A^ Al~^l~ 



where z = cos9 and 9 is the angle between the momentum of A;, baryon and /~ in the dilepton 
rest frame; A(o, b, c) = + b'^ + — 2ab — 2ac — 26c. Following Refs. 2j, |9l|, the differential and 



normalized forward-backward asymmetries for the semi-leptonic decay A;, — > AZ+Z can be defined 



as 



dq- J,^'^q^-Jj'^^- 
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and 

Jo dg^dz + J-1 dg^dz 

Making use of the decay amplitude in Eq. ([5]), the differential forward-backward asymmetry for 
decays of A;, — > A + l^l" can be calculated as 



= 256m3 ^5 A(m^2,m^,g ), (85) 

^ A;, 



If, 

with 



RpBiq^) = 2[{msmA + mfomAj/l - ms{m\ - + q^)f292 + ijUsTriK - mi,mK^^)q^ gl]Yle{Cj^ ^ CIq 
+ [(/2 - g2mK? - glmlM^'^<ClffCl^), 

where we have retained masses for both the legton and strange quark. In the limit of — > 0, 
our results will be the same as that in Ref. [22]. The differential forward-backward asymmetry 
for semi-leptonic decay of Af, — > A/+/^ only depends on the following two combinations of Wilson 
coefficients Re(C7^^^C7*o) and ReiCl^^Cl^), since only terms involving Tr[L^L^] in the differential 
decay width can give rise to one power of "cos0" [22], where and represent the vector and 
axial- vector currents of leptonic sector. The zero position Iq of forward-backward asymmetry is 
determined by the relation 

Re{ClffCl,) = J-^^^^^ f'-'f2 ^ Re(C-//CIo), (87) 

to (/2 - 'mAg2r - iTT'A.ai 

where the small number of strange quark mass is neglected. It is easy to observe that Iq only rely 

on the Wilson coefficient as well as the ratio of form factors. Utilizing the form factors computed 

in this work, we can derive the zero position of forward-backward asymmetry at to = 6.0GeV^ 

without long-distance contributions, which indicates that the forward-backward asymmetry for 

Afe — > A -|- r+T^ could not be zero apart from the resonance regions and end points. However, the 

number of to can shift to 3.5GeV^ using the form factors calculated in QCDSR [34], due to the 

qui t e di fferent predictions on the ratio of from factors in these two approaches. It is mentioned in 



lOOl ] that the zero-position of forward-backward asymmetry is not sensitive to the form factors 
in the large energy limit. However, power corrections contributed by the higher twist distribution 
amplitudes of A baryon can bring about significant effects on the ratio of form factors, which can 
even overwhelm the sign of form factor breaking the heavy quark symmetry. The position of the 
zero of the forward-backward can shift observably for the ratio of form factors with different sign. 
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The distribution of normalized forward-backward asymmetry as a function of with and without 
long-distance contributions is presented in Fig. [5] and Fig. [6l 



O.IS 




a 5 10 15 20 

(a) 




14 16 IS 20 



FIG. 5: Normalized forward-backward asymmetry Aps^q^) as a function of with long-distance contribu- 
tions: (a) for muon, and (b) for tauon. 

It is also useful to introduce the integrated forward-backward asymmetry (Aps) in order to 
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FIG. 6: Normalized forward-backward asymmetry AFB{<f) as a function of without long-distance con- 
tributions: (a) for muon, and (b) for tauon. 
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characterize the typical value of forward-backward asymmetry 2^] 



{Afb) = / AFB{t)dt, (88) 

with t = /m\^, tmin = 4m^/m^^ and tmax = {n^\^^ — ''TT-jJ/^A^- The numerical results of inte- 
grated forward-backward asymmetry are grouped in Table IVIIl where the evaluations in QCDSR 
are also given for a comparison. As can be observed, the predictions on the {Aps) in LCSR are 
typically one order smaller than that in QCDSR due to the quite different predictions on the ratio 
of form factors. The large discrepancy may be recoiled by including the higher conformal spin 
contributions in the sum rules of form factors on the light-cone as well as radiative corrections for 
both two types of sum rules. 

TABLE VII: Averaged forward-backward asymmetry (%) for A;, A + l^l^ with and without long-distance 
contributions based on the form factors calculated in this work using light-cone sum rules and that from 
QCD sum rules 



model 


Ab 


^ A + fi+fi- 


Afe ^ A + T+T" 


LCSR 
(this work) 


without LD 
with LD 


1 99+1.42 
^•^^-0.73 


-0.21 
'-'•"^-0.21 


QCDSR [22] 


without LD 




_3 qq+0.01 



VI. SUMMARY 



The study on rare decay of A^ — > A -|- 7 and A^ ^ A+l^l can serve as the baryonic coun terp arts 



of analysis on i3 — > K*^ 



11 



12 



13 



A-t-7 ana Ah ^ il 

14, i3,Q,[i3,Q, 



19 



2C 



m,B 



[BE 

6d] and A, 



10(] and 

^ AW 



K*rr [51, 

could also be the extension of investigations on heavy baryon decays Af, pW 
64( 1 . Such decays play the role as a corner stone [121] to explore the quark-flavor structure of the 
SM as well as determine its fundamental parameters such as the CKM matrix. 

Although we have achieved inspiring progress in the heavy meson decays with the help of 
heavy quark expansions and factorization techniques, serious and systemical studies on heavy 
baryon decays are comparatively behind on account of the complexity of inner structures for baryon 
systems. In particular, the predictions on the physical observables associating with heavy baryon 
decays can vary even by the orders of magnitude employing different theoretical tools. In this 
work, we explore the soft contributions to the form factors responsible for the A^ ^ A transition 
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in terms of LCSR approach to understand the tremendous discrepancy on predictions of decay 
rate for Af, ^ A + 7 between theoretical methods. More importantly, power corrections from the 
higher twist distribution amplitudes of A baryon to the transition form factors are also investigated 
to the leading conformal spin in detail. It is observed that the higher twist LCDAs almost have 
no influences on the transition form factors reserving the heavy quark spin symmetry, while such 
corrections can result in significant impacts on the from factors breaking the heavy quark spin 
symmetry. In addition, we also find that various models for the distribution amplitudes of A 
baryon can lead to quite different predictions on the transition form factors, which will be tested 
by the experiments on Tevatron and LHC in the future. 

We confirm the conclusion of Ref. 3J, |83] that the A polarization asymmetry of A;, — > A7 
only relies on the relative strength of left- and right- handed couplings between quarks and is 
free of the pollution from the strong interaction. We also discuss some interesting observables in 
phenomenology, such as decay rate, dilepton distribution, polarization asymmetry and forward- 
backward asymmetry of A^ A/+/~. We use the pole model to extrapolate the results of form 
factors calculated in LCSR to the whole physical region in view of the failure of light-cone expansion 
for correlation function in the large momentum transfer region. Our results indicate that the decay 
rate for A^ A/x^/i" is about five times larger than that for A;, A7 due to the long-distance 
contributions from the charmonium resonance region. The polarization asymmetry of A;, — > Afj,^fi~ 
is much sensitive to the long-distance contributions than that of A^ — > Ar+r^, since leading 
resonance J/V' does not contribute to the case of final state with tauon pair. As for the integrated 
forward-backward asymmetry, our results for the magnitude of both muon and tauon cases are 
typically one order smaller than that given in QCDSR due to the quite different predictions on 
the ratio of transition form factors, which can also shift the zero position of forward-backward 
asymmetry from Iq = 6.0GeV^ in LCSR to 3.5GeV^ predicted by QCDSR approach. 
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APPENDIX A: QCD SUM RULES FOR /a, /a, AND Ai 



The sum rules of /a and Ai have been derived in 6J], we would like to collect them in this 
appendix for the completeness of the paper. 

(4^)Vle-^'/""s = I r%(i_;,)5e-/Af|^^_^ r\(i_^)(i_2^)e-/M|^, (AI) 
5 J^2 6 J^2 s 

A{27:Y\lM\-^''^'l = 1 r [(1 _ x^){l - 8x + x^) - I2x^ In x] e-'/^^Bds 

+ ^ r\l _ xfe-^/^^Bds - - a^e-^'/^B, (A2) 

4 J ml 3 

where x = ml/ s and is the mass of strange quark. The values of the non-perturbative conden- 
sate at scale ^ = IGeV are given by 

a = -{2TTf{qq) = 0.55 GeV^ 

b = (27r)2(^G2) = 0.47 GeV^ . (A3) 
vr 



The sum rules for can be read from Eq. (|Aip . only with the following replacement: 

/a /At, ms rrib, s% s^^. (A4) 

With the threshold value = 1.6GeV^, we can arrive at the numbers of /a and Ai as 

= 6.0+°:^ X lO-^GeV^, Ai = -1.31°;^ x lO'^GeV^ (A5) 

within the Borel window M| G [1.0, 2.0]GeV^. In the same way, the result of /a^ can be derived 
as /Afc = 3.9l^Q2 X 10^'^GeV^ with the selected threshold parameter s^^ = 39 it IGeV^ and Borel 
mass M| G [2.0, 3.5]GeV^ 
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